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Abstract 

We investigate the infrared (IR) effects of Lorentz violating terms in the gravitational 
sector using functional renormalization group methods similar to Renter and collaborators 
[T]. The model we consider consists of pure quantum gravity coupled to a preferred 
foliation, described effectively via a scalar field with non-standard dynamics. We find 
that vanishing Lorentz violation is a UV attractive fixed-point of this model in the local 
potential approximation. Since larger truncations may lead to differing results, we study 
as a first example effects of additional matter fields on the RG running of the Lorentz 
violating term and provide a general argument why they are small. 

1 Introduction and motivation 

One of the major unsolved problems of theoretical physics is the construction of a physically 
acceptable and predictive UV completion of quantum gravity. The perturbative nonrenor- 
malizability of gravity hints that its quantization requires new physics, such as a preferred 
foliation modeled, e.g., in Einstein-aether theory |2j. Preferred foliations gained a lot of atten- 
tion recently after Hofava's proposal [3] of a perturbatively renormalizable quantum theory 
of gravity. Hofava's model is a higher derivative gravity theory that avoids troublesome 
unitary ghosts by retaining second derivatives in time while allowing higher derivatives in 
space. The preferred foliation can be encoded [5] in a Stiickelberg-type field (j) which defines 
an irrotational unit time- like vector field = 4',fi/ \/9'^''4',u<P,p, where the preferred folia- 
tion is defined by constant cf) surfaces. The extrinsic curvature of the preferred foliation is 
K^u = ^niSfMu + n^n^), where Cn denotes the Lie derivative in direction of n. The low energy 
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limit of Hofava's theory contains in addition to the usual Einstein-Hilbert action also the 
Lorentz violating action given by 
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^\fW\ + higher derivative terms , 
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where Gn denotes Newton's constant and K = g^^^K^^,. 6 is a new coupling constant that 
vanishes in general relativity. A possible K^'^K^^-term is absorbed in Hoi^ava's definition of 
the speed of light. The modification ([1]) is a low energy effective field theory, describing the 
effects of quantum gravity at a scale k where higher derivative terms are negligible. In order 
to bring this theory into contact with observation we use the renormalization group equation 
to evolve this action to IR scales at which experimental tests for Lorentz violation are being 
performed. 

A Wilsonian approach to the renormalization group equation for ([1]) leads to considering 
nonperturbative renormalization of Einstein-aether theories, which is very involved. A much 
simpler on-shell equivalent to low energy Hof ava theory has been found in [5] in terms of the 
cuscuton field: 



This model generates the term dynamically and, given sufficient boundary conditions, 
leads to the dynamical emergence of a constant mean curvature foliation. The quadratic 
cuscuton coupling oi (which has dimension mass~^ in our normalization) is related to the 
Lorentz violating parameter b by 



The observational constraints on the gravitational Lorentz violating parameter b are much 
weaker than the bounds on Lorentz violating parameters in extensions of the standard model. 
A cosmological bound constraining —b = 0.003 it 0.014 at 95% CL can be found in [6]. 
Notice that the Lorentz violating parameter b is dimensionless, so there is no power counting 
argument for its running. 

The advantages of Renter's programme, which achieves to be at the same time nonpertur- 
bative and background independent, suggest to use this Wilsonian framework for the inves- 
tigation of Lorentz violating models. As a start into the Wilsonian investigation of Lorentz 
violating models of this type, we will consider an adapted version of the cuscuton model ^ 
generalized to an arbitrary even power-series potential V{(p) and discuss the /3-function of the 
quadratic cuscuton coupling in presence of pure quantum gravity. The main result is that 
in the local potential approximation the /3-function of the Lorentz violating parameter b has 
a UV attractive Gaussian fixed-point 6 = 0, i.e. the /3-function of the absolute value |5| is 
negative in the vicinity of this fixed-point. Fortunately, the corresponding IR growth of b is 
bounded due to an IR decreasing /3-function /3(6), if we assume a sufficiently semi-classical 
RG fiow in the gravitational sector. We show that coupling additional matter fields to this 
theory has only minor effects on the RG flow of the Lorentz violating parameter, at least for 
small matter couplings and Lorentz violations b. Notice that since the cuscuton model differs 
from Hofava's theory, both through off-shell terms and higher derivative contributions, one 
can not apply these results to Hofava's theory itself, but one should consider them as a toy 
model to gain a first insight. 
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The outline of this paper is as follows: In section [2] we provide the background on Wet- 
terich's functional renormalization group equation necessary in order to be self-contained. 
Section [3] explains the basic idea of IR attractive gauge symmetries using a toy model. In sec- 
tion m we calculate explicitly the flow of the parameter b within the cuscuton model coupled 
to quantized gravity. The introduction of matter to this model is discussed in section [5l We 
give an outlook towards further investigations in this direction in section [6l 



2 The functional renormalization group equation (FRGE) 

In this section we present a brief review of Wetterich's functional renormalization group 
equation (FRGE) [7j to provide the prerequisites for this paper using a Euclidean scalar field 
theory; detailed introductions can e.g. be found in [8] [9] [10] [11] and reviews on the application 
to gravity in [T^Pj[14j. 

The starting point for formulating the FRGE is the partition function defined by the 
regularized Euclidean path-integral 

exp{Wk[J]) = J Vax e-^W-A5,[x]+J.x ^ (4) 

where A denotes an overall UV cutoff, S the bare action (at the cut-off) and ASk an IR 
regulator that suppresses momentum modes with < k"^ by giving them a mass of order k 
while effectively vanishing for modes with > k"^ . A typical suppression term is given by 

A5,[x] = \jd'^ x{x)Rk{-U)x{x) = ^ / x{p)Rk{p')x{-p) . (5) 

The effective average action is defined via the Legendre transformation and subsequent 
subtraction of the IR regulator 

' d^x J(x)</)(x) - Wk [J] - ASk M , (6) 

where (/>(x) = {xi^)) = 5Wk[J]/6J{x) is the expectation value of the field x- To calculate the 
scale derivative dt ■= kd^ of F^, we take the scale derivative of ([H]), substitute with dU 
and reexpress the RHS in terms of variations of F^, resulting in Wetterich's FRGE 



Rk-{Tf\ct>]+Rk) ' , (7) 



(2) 

where dot denotes dt and F^ denotes the second variation of F^, w.r.t. the fields. The trace 
extends over all fields. 

The usual effective action F is attained in the limit /c — )• 0, since AS^ vanishes in this limit. 
The usual bare action is attained in the limit A;, A — )• 00, because ASk diverges in this limit, 
making the saddle point approximation for ^ exact. The effective average action F^ is thus 
an interpolation between the usual bare action S (for A — )• 00) and the usual effective action F 
and effectively describes physics at the scale k, because its tree-level accurately describes the 
effects of momentum modes above k. This allows one to derive an effective action without 
explicitly solving a path-integral by specifying an initial condition Tk=A and using ([Tj) to 
evolve Ffc to the physically relevant scale k. 
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To calculate /3-functions from ([7]), one expands in terms of field monomials in the flow 
equation and equates the coefficients of the field monomials on the LHS with the corresponding 
coefficient on the RHS. However, the trace on the RHS of the flow equation produces in 
general an infinite number of field monomials, meaning that one can not stick with an action 
of a particular form, but one is forced to make the most general ansatz compatible with 
field content and symmetries. Using the most general ansatz for the effective average action 
however yields an insurmountable calculation, so for practical purposes one will have to make 
an ansatz by selecting the "most important" monomials and expand the trace on the RHS 
only within this ansatz. The strength of the flow equation is that one can use expansions that 
are not necessarily perturbative, in particular an expansion in the number of derivatives of 
the field operator. The lowest order of the derivative expansion is the local potential ansatz, 
which reads for a standard scalar field 

.4. /I 



d^x(^-a^09'^</> + l^fc((^)j . (8) 

Practically, one will consider only the lowest dimensional monomials, so imposing a Z2 sym- 
metry, we may consider 

^.(0) = ^^^ + ^^^ + ... • (9) 

The focus of this paper is the renormalization of the cuscuton field coupled to gravity, which 
we will do in the local potential approximation. For this we have to give a generalization of 
the FRGE methods so gravity can be included. An important ingredient are background-field 
methods, as used in various applications to gauge and gravity theories. This is particularly 
important in gravity, since coarse graining in a background independent theory needs a split 
of the metric g^^ = g^u + h^^, into an arbitrary background metric gfj,^, and fiuctuations h^^, 
which are not assumed to be small. Following Renter [1] , we use the background Laplacian to 
define the momentum scale k of the fluctuation field. In contrast to the perturbative graviton 
expansion however, one does neither assume a fixed background nor that the fiuctuations are 
small; it is hence not necessary to make an expansion in terms of gravitons for the price that 
the effective action of gravity depends on both, g^^, and /i^^. 

The trace on the RHS of the flow equation can then be expanded background indepen- 
dently in terms of curvature invariants using the heat kernel expansion. The first two terms 
in this expansion yield the invariants occurring in the Einstein-Hilbert action. The Einstein- 
Hilbert truncation with harmonic background field gauge fixing reads 

TEu[g, h] = j d^x^[-R[g] + 2A) 

+ j d^x^DP (^h^p - r^D^ - + Shosts , (10) 

where h = g^'^h^^ and = 327rGAr. We did not give the explicit form of the ghost action, 
since it does not explicitly contribute to the running of the quadratic cuscuton coupling in 
our truncation. 

The Rk and terms appearing in the FRGE ([7]) include only the functional derivatives 
with respect to the fiuctuation fields. The background fields set the scale k present in the 
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regulator function, i.e. we define for h^y, as well as for all other dynamical fields, the regulator 
action with respect to g^u and the corresponding Laplacian □ 

A5fc[^, h] = K^ j d^x^ V ^^"(5, -□) V • (11) 

The background gauge condition is satisfied for vanishing fluctuations, so we may, after the 
FRGE has been derived, set the averaged fluctuations to zero, i.e. hfj^iy = 0. Then the back- 
ground field g^u can be interpreted as the averaged quantum metric field. 

A useful trick [15] in the evaluation of the trace on the RHS of the flow equation is to 
obtain a formal expression iTr(0[(/)]) = so the flow equation reads AjOi['/'] = ^[4>]^ 

where Aj denotes a coupling and Oi[(j)] a field monomial. It is now useful to insert families of 
field configurations that (1.) project onto the truncation on the LHS, (2.) identify the terms 
in the truncation uniquely and (3.) simplify the evaluation of F, e.g. constant matter fields 
to project onto a local potential ansatz or flat metrics to project onto monomials that do not 
involve curvature invariants. 



3 IR symmetries: General scenario 

Let us assume that we integrate out short distance physics above a very high scale A and 
hence produce an effective average action rA[</>j] with field content {(j)i}iex. Let us furthermore 
consider a gauge group Q acting on the field content so the gauge invariant field content 
is {(pi}/Q. The question that we want to consider is whether this particular gauge symmetry 
is IR attractive or whether one needs fine-tuning to retain gauge symmetry in the IR. To 
answer this question we use the exact renormalization group equation and evolve Fa from 
A down to an effective scale k and observe whether or not the /3-functions of the couplings 
between physical and gauge degrees of freedom are positive for a finite range of couplings. 
Since the exact renormalization group for the effective average action requires us to consider 
the most general action that is compatible with our field content and symmetries, we have to 
consider a truncation that contains all couplings that are "low dimensional". Similar ideas 
have been investigated in [16j [17J [18j . 

Let us consider an ultralocal toy example consisting of a real scalar multiplet {(pi}'f^i with 
an ultralocal gauge group acting as 

^i{x) ^ X{x) + (j)i{x) . (12) 

Let us change to the variables 

(fi = 4>i+i for 1 < i < n , , , 

where {(pi}^Zi physical and ^ is pure gauge. To break this gauge invariance in a controlled 
manner let us introduce an additional "Stiickelberg field" x transforming as x(^) n\{x) + 
xix) under gauge transformations, so the most general local potential ansatz with standard 
kinetic term allowed by gauge symmetry is 

Lfc = / d^x \d^v,d^^, + ]^d^{^ - x)d^{^ -X) + Vk{^^, f - x) ) • (14) 
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Ultralocality allows us to fix a gauge x = 0, making $ physical in this gauge such that the 
original gauge symmetry means that $ decouples from the theory. Let us simplify the dis- 
cussion by considering n = 2 and a truncation of the local potential to at most dimensionless 
couplings and let us impose an individual Z2 symmetry, such that the truncation becomes: 



k nil. 9 k M, ^9 ah A hh Ch ^, 



2 ^ 2 4! 4 4! 

The /3-functions can be extracted from 



(15) 



«M|) = -2Ml - ^ + ^^^) , (16a) 

R(h ^ - ^ ( "^fe^fc I bkCk 

^^''=^"(47r)2 V(l+m2)3 + (l + M2)3 

+ n + '3 ^ (16b) 

^^''^^ = ii^{-ir^^^jr^) • ^''^^ 

We see that the /3-function of the dimensionless mass M| is always negative whereas the 
/3-functions of the four-point couplings and are always positive if ak,bk,Ck > 0, meaning 
that the decoupling limit Mk large and bk — )• is IR attractive in this regime. This implies 
that the renormalization group flow drives towards original gauge symmetry, so no fine-tuning 
is necessary. 

The opposite scenario is more familiar: If we consider a generic gauge theory, but forget 
about gauge symmetry, then there are in general many low dimensional field monomials that 
are forbidden in the gauge theory that may enter the theory without gauge symmetry. It is 
very unlikely that the unphysical degrees of freedom will decouple from the physical degrees 
of freedom, since this requires that all new field monomials are RG irrelevant, whereas low 
dimensional field monomials have the tendency to be RG relevant. 

The most robust predictions of Wilsonian renormalization come from the discussion of 
fixed-points, which corresponds to the determination of universality classes. The IR-symmetry 
scenario can be seen at a glance when one considers the Gaussian fixed-point at which 
the masses are relevant and the four point couplings are (marginally) irrelevant. The IR- 
symmetries in the vicinity of any fixed-point are given by the symmetries of the relevant, 
marginally relevant and completely marginal directions at the fixed-point. Notice that i.g. the 
IR-symmetry may be affected by a nontrivial wave function renormalization. However, this 
is not the case if we are dealing with a Gaussian fixed-point due to the following argument: 
If we introduce wave function renormalizations Zi , Z2 in equation (|15p then the /3- function 
of the symmetry breaking coupling B = Z\Z'2b is /3(i?) = {r]\ + r]2)B + ZiZ2(3{b), where 
r]i = dtZi. We obtain that B has a Gaussian fixed-point if and only if b has one and that 
dsPiB) = db(3{b) at the Gaussian fixed-point b = B = 0. Note that the stability of a non- 
Gaussian fixed-point may be affected by r/j. The running of Zi was however minimized in 
test models by using the optimized cut-off [T9], so to gain a first insight, one may neglect the 
effect of r]i, even for the discussion of generic fixed-points. 
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4 Explicit renormalization of the cuscuton model 



In this section we discuss the renormahzation of an adapted version of the cuscuton model 
[SI [6] coupled to the Einstein-Hilbert truncation of general relativity in order to investigate 
the quantum behavior of the term in the Lorentz violating action ([U). We make a general 
Z2-symmetric power-series potential ansatz for the effective action 

Tk[g,K<P]=T^n[9M + I d^x^ |^^5'^-</>,^</.,. + + E (5)! '^'") ' ^^^^ 

including the quadratic cuscuton term ([2]) for n = 1. Note that, in order to have a background 
independent FRGE, we introduced A = fc^A with an arbitrarily small A in the kinetic term. 
This parameter is treated as an external parameter and defines a family of quantum field 
theories. We show that the limit A — )• is well defined after the FRGE has been derived, 
leading to the cuscuton model. The sole purpose of this parameter is to regularize the cuscuton 
propagator without introducing a cuscuton background, as it will become clear below. We 
will treat gravity within the Einstein-Hilbert truncation and focus on the renormalization of 
the coupling constants k?, A and an- 

We start by collecting all terms required to evaluate the RHS of the fiow equation ([7]). 
The second variations T^ are given in the appendix [Bj Without loss of generality, we may 
insert families of fields that project onto our truncation and which allow us to distinguish 
the field monomials therein. For extracting the /3-functions of the coupling constants a„ it 
is therefore sufficient to insert a flat Euclidean background metric Qf^uix) = 6^i, as well as a 
constant cuscuton field 0(x) = (p. Note that this is only a convenient choice that simplifies 
our calculations; the /3-functions constructed therewith are generally valid. 

As regulators {Rk)AB for the individual fields A, B £ {0, /i/^i/} we use the properly nor- 
malized optimized regulator [20] and define 



{Rkh<t, = \{D + k^)e(l + ^) , (18a) 

{RkThl'"' = 2k^K^^'"^{U + A;2)G(1 + ^) , (18b) 
where all indices are raised and lowered with the background metric and where 

The regularized second variations (F^ -|- Rk)AB can be split into a (j) dependent part 
{Fk)AB and cf) independent part {Gk)AB- The inverse is then given by the geometrical series 



(f(') + Rk)-j,], = {Pk)AC ■ £(-Ffc • PkYcB , (20) 



^(2) 

-r ^i-kJAB — V-f k)AC ■ 

1=0 

where {Pk)AB '■= (G'fc)^]^ are the regularized propagators given by 



Ajy-p.i'pu 



2k2(A;2-2A 



A 

2k^ 
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For ^0 and k^-2A^0 we can always choose A sufficiently small in order to expand (j2ip 
in A. Thus the propagators are power-series in A with non-negative powers. 

It turns out that the cuscuton does not contribute to the flow of the coupling constants 
an in the limit A — >• 0, since all internal cuscuton propagators come with A^ and all field terms 
{Fk)AB come with non-negative powers of A as well. Even traces including {Rk)4,<f> ~ A^^ 
do not contribute to the flow of the cuscuton potential, since they occur in the following 
combination 



Tr 



{Rk 



(22) 



The cuscuton only contributes to the flow of the gravitational sector, i.e. k and A, via the 
trace 



1 



Tr[{RkU^-{Pk)^^] =0(AO) . 



(23) 



In the limit A — )• this contribution to the renormalization of n and A is equal to the 
contribution of a canonical massless scalar fleld and does not depend on the cuscuton potential. 
Thus we flnd that in the limit A — >• only gravity contributes to the /3-functions of a„. In 
particular, only one term is contributing to the /3-function of the quadratic term ai 



with the fleld term 



-It, 

2 



■klhh ■ {Pk)hh 



k hh 



k hh 



(24) 



{F, 



k hh 



(25) 



This leads to the dimensionless /3-function for ai = k'^ai: 

/ 

/3(ai) = ai 2 + 



/3(k^) + 8k 



2^2 



\ 



(26) 



where we introduced the dimensionless coupling constants Aj deflned by Aj = /c'^'^^'Aj. The 
Lorentz violation is governed by the parameter b = (4aiK^)~^, see ([3]). Its /3-function is 



m 



+ 



1927r2K4 



1 - 2A 



(27) 



The /^-functions of and A are equal to Renter's /3-functions [Ij plus the additional contri- 
bution of a massless free scalar due to ([23]) . which can e.g. be found in [13] [2T]. 

We observe that the /3-function of h has a Gaussian flxed-point. Assuming real flxed-point 
values K=|,,A* we see that dh(3(b) < 0, meaning that the (6 = 0, k = k*, A = A*) flxed-point 
is UV attractive. Assuming a not necessarily Gaussian flxed-point k*,A=k, there appears to 
be the possibility for (3(b) to vanish for arbitrary values of b due to vanishing of the bracket 
in (j27p . This bracket can however not vanish for real values of , A^^ , meaning that the only 
physically acceptable flxed-point for bisb = 0. This means that the scenario of IR-attractivity 
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of Lorentz symmetry, as explained in section [3l is not realized in the cuscuton model. This 
does however not mean that this scenario may not be realized in models with a different 
implementation of the Lorentz violation. 

If we make the simplifying assumption that the /3-functions of and A behave classical 
in some finite interval of the scale k, meaning they have no anomalous scaling 

/3(k2) = -2^2 , /3(A) = -2A , (28) 

then we obtain for (j27p 

/3(6) = , , (29) 

where t = Ink and denote the "initial conditions" at t = 0. The physics described 

by this /3-function (j29p is as follows: The absolute value of the Lorentz violating parameter 
grows in the IR, i.e. as t — )■ — oo. This increase is bounded due to an exponentially decreasing 
/3-function. Thus one can in principle force b to experimentally valid values by a suitable 
fine-tuning, but the phenomenologically optimal scenario of a highly Lorentz violating UV 
theory which flows down to a sufficiently Lorentz invariant IR theory can not be realized 
without seemingly unnatural assumptions or an unexpectedly large wave function renormal- 
ization, which we find unlikely due to our use of the optimized cut-off [19]. The practical 
calculation of such a wave function renormalization would however require new techniques for 
the evaluation of the traces, since the cuscuton kinetic term can not be generated with the 
derivative expansion nor with a vertex (number of fields) expansion. 



5 Adding matter fields 

Since the renormalization of the quadratic cuscuton term coupled to pure gravity obtained in 
the previous section does not lead to IR attractivity of Lorentz symmetry, we will now study 
whether this can be achieved in presence of additional matter fields. The simplest standard 
matter model is an additional scalar field x with canonical kinetic and mass term, together 
with a coupling to the cuscuton field via the lowest-order interaction 

.^int = 5mat j d'^X^X^^''^ ■ (30) 

Note that even if this interaction would not be present at some scale, it would be induced 
radiatively via the graviton interaction at a different scale. It is thus natural to consider such 
an operator. We find that one additional trace contributes to the flow of b, namely 



2 



iRk)xx ■ (^fc)xx ■ i^k)xx ■ (^fc)xxj ■ (31) 
The terms involved in this trace are given by 

{Rk)^^ = {a + k'mi + ^) , (32a) 

(P^O.v = 1^ ■ (32b) 

A 

(Ffc)^x = 2(7mat0^ . (32c) 
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and the resulting /3-function for b is given by 



k2 



5 P{k'^) + 8k2 




= -b 



+ 



27r2(l + m2)2 



(33) 



V 



We observe that this /3-function can in principle be positive for a suitable choice of b and 
5mat- Since we are phenomenologically interested in small Lorentz violation |6| ^ 1 and 
small matter couplings |gmat| ^ 1 (the interpretation is invalid at the scale at which the 
dimensionless matter couplings are not small), we obtain that matter i.g. does not change the 
sign of the /3-function in our parameter range. On the other hand ()33p leads to a condition 
for a non-Gaussian fixed-point 6 = 6* 



This means that if one had found a gravity-matter system with non-Gaussian fixed-point, then 
there could be a non-Gaussian fixed-point for the gravity-matter-cuscuton system. However, 
even if one found such a system, one would still generically have a negative /3-function in the 
phenomenologically required regime, where |6| <^ 1 as well as l^matl 1- 

Similar effects can be obtained by introducing other matter fields, if they give rise to 
positive contributions to the /3-function of 6. In particular, all matter effects will enter with 
^2 and will be required to couple weakly to the cuscuton field, so the interpretation holds. 

6 Conclusions and outlook 

Motivated by the possibility of constructing UV complete theories of quantum gravity in 
explicitly Lorentz violating theories with preferred foliation, we investigated the low energy 
effect of a preferred foliation in pure quantum gravity. The particular model that we investi- 
gated was a version of cuscuton theory that has recently been shown to be on-shell equivalent 
to Hofava-theory at low energies. We used a truncation ansatz consisting of Einstein-Hilbert 
gravity coupled to a cuscuton field with local Z2-symmetric potential in Wetterich's exact 
renormalization group equation and applied Renter's background field methods for gravity. 
We were particularly interested in the nonperturbative and background independent renor- 
malization of the quadratic cuscuton term which is related to the Lorentz violation. The 
nondynamical nature of the cuscuton field leads to very stable predictions for the renormal- 
ization of the cuscuton potential, which can be traced to the fact that only gravity propagates. 
We find in particular: 

1. The /3-function of the quadratic cuscuton coupling is universal in the sense that it is 
independent of the rest of the cuscuton potential parameters. 

2. The only fixed-point of the Lorentz violating parameter is a UV-attractive Gaussian 



3. The inclusion of generic scalar matter has very weak effect on the UV-attr activity of 
the Gaussian fixed-point. 




(34) 



fixed-point. 
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4. Although having a quadratic term, the effect of the cuscuton on the renormahzation of 
Newton's constant and the cosmological constant is equivalent to the effect of a massless 
scalar field. 

The physical interpretation of these results is that there might be some Lorentz violation at 
very low energy scales, such as the Hubble scale, that may lead to observable cosmological 
effects. But using current cosmological bounds on Lorentz violation at the Hubble scale we 
expect no observable Lorentz violation in the high energy regime. This result contradicts the 
optimal scenario in which the renormahzation group drives an initially large Lorentz violation 
in the UV theory to an unobservable value in the IR, hence if Lorentz violation is necessary to 
the UV theory one needs (within the universality class of the cuscuton model) an additional 
mechanism to tame this Lorentz violation in the IR or an unexpectedly large wave function 
renormahzation. 

To gain further insight into the effect of preferred foliations one has to investigate whether 
different gravity models with preferred foliations lie in the same universality class as the 
cuscuton field. Within the cuscuton theory itself one could gain further insight by performing 
calculations in truncations beyond the local potential ansatz, in particular the wave function 
renormahzation. The investigation of the cuscuton wave function renormahzation can not be 
calculated with standard techniques, because it is neither accessible in a derivative expansion 
nor in a vertex expansion. 
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A Toy model calculation 

(2) 

Using the same regulator for all fields, the fluctuation matrix + -Rfe of (|15|) can be 
split into F + G with: 

/ R,+p^ + k^ml A 



To project the RHS of Wetterich's equation onto the local potential ansatz, we may insert 
constant fields, and extract the /3- functions for the masses from the trace —^Tr (PRPF) and 
for the four-point couplings from ^Tv{PRPFPF) respectively, where P = G^^ . Using the 
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2 

optimized cut-off Rk = {k"^ — p'^)@{l — we find for constant fields: 

-iTriPRPF) = -6&/d^x((^+(^)^^ + ((^+(^)cI> 

^LT^r-nr-r r-r ) - ] a xy 2(i+m^)3 ^ ^^x+miY{\+Ml) ^ (i+m^)(i+Mfc^)2 + 2{i+M,^)3 

(36) 

B Second variation of the cuscuton model 

We give tlie complete second variation of the effective action (fT7|) . It is obtained by plugging 
in the background field expansion g^^y = gnu+eh^^ and cp = (p+ef into the effective action (fT7|) 
and collecting the term^. For a better readability we give independently the individual 
parts of the second variation. 

The second variation of the kinetic cuscuton action is given by 
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The second variation of the cuscuton potential part is 

(2n)! ' 2 (2n - 1)! ' 2 (2n- 2)! 



pot ^ " I \ S zl / ('9'nV 9 ('9ti - m 9 ('9r7 - 9V I ^ ^ 



n=l 



The second variation of the Einstein-Hilbert part can be found in \T\ for general background 
metric fields. Insertion of a flat background metric into the general variation yields: 

rgL = J d^x h^.K^^^'P^U + 2A) V , (39) 



where K^'"/'^ was defined in (|T9|) . 

The functional derivatives with respect to the individual fluctuation fields {/, h^^} of the 
sum of the expressions above enter the RHS of the FRGE ([7]). 
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